Novel Extrapolation Method in the Monte Carlo Shell Model 
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We propose an extrapolation method utilizing energy variance in the Monte Carlo shell model 
in order to estimate the energy eigenvalue and observables accurately. We derive a formula for 
the energy variance with deformed Slater determinants, which enables us to calculate the energy 
variance efficiently. The feasibility of the method is demonstrated for the full p/-shell calculation of 
56 Ni, and the applicability of the method to a system beyond current limit of exact diagonalization 
is shown for the p/+<?9/2-shell calculation of 64 Ge. 

PACS numbers: 21.60.Cs, 27.50.+e, 24.10.Cn 



The shell model (SM) calculation has been very suc- 
cessful in understanding the nuclear structure on the ba- 
sis of nucleons interacting via the nuclear force. The 
conventional, standard solver for SM calculations is the 
exact diagonalization of Hamiltonian matrix in a given 
model space. Recently, the SM calculation plays an in- 
dispensable role especially in studying neutron-rich ex- 
otic nuclei, including beta-decay properties on r-process 
nuclei {e.g. [JIH). For such studies, the model space of 
the SM calculation should contain some intruder orbits 
in addition to one major shell. In this case, the dimen- 
sion of its Hilbert space is often explosively large and the 
practical calculation is infeasible. Overcoming such a dif- 
ficulty is a crucial challenge for modern SM calculations, 
where much effort has already been directed {e.g. 043 )■ 
The Monte Carlo shell model (MCSM) Q is one of the 
methods which aim at surpassing the limit of the conven- 
tional diagonalization Q and have succeeded in realistic 
applications. 

The MCSM has been formulated by combining 
auxiliary-field quantum Monte Carlo and diagonalization 
methods [9(. The MCSM yields the resulting wave func- 
tion as a linear combination of a relatively small number 
of deformed-basis wave functions. While the convergence 
pattern of the energy eigenvalue as a function of the basis 
number suggests the validity of the approximation, the 
convergence is, in many cases, not fast enough to estimate 
the exact energies accurately. This is a long-standing 
problem in the MCSM. The same problem also occurs in 
the conventional SM calculations when the model space 
is truncated. 

In the case of the conventional SM calculations with 
truncation, the approximated eigenvalue seems to de- 
crease exponentially as a function of the basis number. 
As an empirical trial, the exact energy can be guessed 
by an exponential extrapolation 10], though this tech- 
nique cannot be applied directly to the MCSM. In this 



paper, to estimate the exact energy eigenvalue, we con- 
sider another novel method free from such convergence 
patterns. 

Recently an extrapolation method utilizing energy 
variance to estimate exact energy eigenvalue has been de- 
veloped [ll| ■ Because this method is expected to be valid 
independently of the representation of the basis function, 
its application to the SM is of interest. In spite of efforts 
for such applications [IH, [l3j], its full-scale application 
has been infeasible due to the limitation of computer re- 
sources. In the present work, by deriving a new formula 
for the expectation value of the Hamiltonian squared, 
such an extrapolation is made feasible. 

First, we briefly review the framework of the MCSM. 
We use a general two-body interaction as: 
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where cj denotes a creation operator of single particle 
state i. In the present work, the MCSM wave function 
is given as a linear combination of angular-momentum- 
projected, parity-projected deformed Slater determinant 
wave functions, 
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where PjJ K is the angular-momentum and parity projec- 
tor, and N is called the MCSM dimension. Each \ip n ) is 
a deformed Slater determinant, 



ku=n e^ 
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where |— } denotes an inert core. The coefficient D t - n " > 
is selected from many (roughly one thousand) candi- 
dates generated stochastically utilizing the auxiliary field 
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Monte Carlo technique. The coefficient ^ is deter- 
mined by the diagonalization of the Hamiltonian matrix 
in the subspace spanned by projected Slater determi- 
nants, P^ K \ip n ). This diagonalization also determines 
the energy En = (^nIHI^n), as a function of N. In 
principle, we increase N until En becomes converged. 

Next, we introduce the energy-variance extrapola- 
tion into the MCSM (MCSM-extrapolation method) fol- 
lowing the idea of Ref. [El]. The MCSM provides 
us with a successive sequence of the wave functions 
^1), 1^2), • • • , \^n), ■ ■ ■ ■ For each N, we evaluate en- 
ergy variance as, 

(AH 2 ) N = (* N \H 2 \* N ) - (V N \H\y N ) 2 , (4) 

and plot the energy En as a function of its variance. 
As we increase N and improve the approximation, the 



resulting energy approaches the exact energy, and the 
corresponding energy variance approaches zero. These 
values are fitted by a second order polynomial, and the 
energy is extrapolated to the limit of (AH 2 ) — > in 
the same manner as other applications of energy- variance 
extrapolation [lS^ . 

The obstacle in the implementation of the MCSM- 
extrapolation method was the large amount of compu- 
tation to evaluate ((f>\H 2 \ip), where \<j>) and \i/j) are de- 
formed Slater determinants. If we regard H 2 as a general 
four-body operator, the evaluation of the matrix element 
consists of the eightfold-loop summation of the 24 terms 
of products of four generalized one-body density matri- 
ces, pij = (<j>\CjCi\ip)/ (<f>\ip). In the present work, thanks 
to the separability of H 2 , the evaluation of the matrix 
element is formulated as: 



^2 V Wl((l - P)ka(l - P)W ~ (! - P)lM - P)k[i) V ^js(p 7 iP5j ~ PSiPjj) 



i<j,a<j3 \k<l 



+Tr((i + T)(1 - p)(t + r)p) + ( Tr(p(t + -r)) 
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with Tik — J2ji v ijkiPij ■ The trivial summations and 
their indices for the matrix products are omitted. The 
first term in Eq.Q is written as a product of two matrices 
as the first term on the right-hand side. This factoriza- 
tion reduces the eightfold loop into a sixfold loop and 
decreases the computation time drastically. 

Now, we apply the MCSM-extrapolation method to 
56 Ni with the p/-shell and the FPD6 interaction [3]. 
The m-scheme dimension of 56 Ni reaches 1.0 x 10 9 . The 
present work was performed using the newly developed 
MCSM code [15l| . which enables us to run it on latest 
supercomputers. 

Figure [1] shows the MCSM results of the ground-state 
(J v = + ) and the first-excited-state (J' r = 2 + ) ener- 
gies as functions of the MCSM dimension. These ener- 
gies show good convergences, but slight differences from 
the exact values remain. We will show how these gaps 
are filled by the extrapolation method later. The en- 
ergy by the current MCSM calculation is -203.161 MeV 
with N — 150, while the past results of the MCSM were 
-203.100 MeV in 1998 [l|, and -203.152 MeV in 2001 
Over a decade, progress in the method and in com- 
putational power has gradually improved the precision of 
the MCSM. Nevertheless, we still find 37 keV error from 
the exact energy, -203.198 MeV. Note that the MCSM 
error of the l\ state is the same order of magnitude. 

Figure [2] shows the En as a function of (AH 2 )n pro- 
vided by the MCSM wave function. We fit the MCSM 
points of E N against (AH 2 ) N with 10 < TV < 150 by 
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FIG. 1: Convergence patterns of the ground and first excited 
states of 56 Ni in the p/-shell. The solid circle and the triangu- 
lar symbols denote the MCSM results of J = + and J = 2+ , 
respectively. The dashed lines show the exact values by the 
diagonalization method. 



quadratic curve, and extrapolate the MCSM results to 
(AH 2 ) -> 0. The extrapolated energy is -203.198 MeV, 
which agrees with the exact one within 1 keV. Here, we 
excluded the first nine points of En for the quadratic fit 
because the extrapolation method assumes that approx- 
imated wave functions are sufficiently close to the true 
eigenstate. Moreover, the MCSM points of N < 10 show 
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FIG. 2: (Color online) Second-order extrapolations of the 
ground-state energy into a zero energy variance of the J n — 
+ ground state of 56 Ni in the p/-shell. The filled symbols, 
open symbols, solid red line, and dotted blue line denote the 
Em of MCSM, the results of the diagonalization method with 
PHT, and their second-order fits, respectively. The exact en- 
ergy is also shown by open symbols on the y-axis. The inset 
shows magnified view around (AH 2 ) ~ 0. 



comparably large fluctuation due to stochastic procedure 
and should have strong dependence on the initial states 
of stochastic sampling. 

For comparison, we also show another extrapolation re- 
sult for the conventional SM calculation with the particle- 
hole truncation (PHT) in Fig. [2l The configuration 
of the PHT is (0/ 7/2 ) 16 -*(0/ 5/2 , 1^3/2,1^1/2)* with t = 
3,4,5,6,7, and the practical calculation was performed 
by the MSHELL code ftfi. These energies and their vari- 
ances are also fitted by a quadratic curve in the same 
manner as Ref. [Hj]. While both the MCSM and the 
PHT calculation succeed in reproducing the exact en- 
ergy well, minor deviation can be seen in the inset of 
Fig. H The extrapolated energy with PHT is -203.217 
MeV, and its discrepancy with the MCSM and the ex- 
act energy is 19 keV. Note that we discuss precision in 
the unit of a few keVs, while previous works using the 
energy- variance extrapolation provided the precision of a 
few tens or a hundred keVs [l2|, EH • 

An advantage of the MCSM for the extrapolation 
method is that the MCSM provides us with the sequence 
of many (more than 50) successive approximate wave 
functions simultaneously. It provides us with good statis- 
tics for the extrapolation. On the other hand, the con- 
ventional PHT scheme yields only 6 points in the case of 
56 Ni, for example. 

In order to test the applicability to larger systems, we 
assume that the MCSM result with TV < 50 is available 
in the ground state of the 56 Ni case. In practical cal- 
culations, the TV is often limited so small that the Em 
cannot reach good convergence. The MCSM result with 
TV = 50 is E N = 50 = -203.115 MeV, which is worse than 
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FIG. 3: (Color online) (a) Second-order extrapolations of 
the energies of J n = Of ,2f ,0^, and 0^ states of 56 Ni in the 
p/-shell. The notation is the same as Fig. f2] (b)First-order 
extrapolation of the occupation number of the O/7/2 orbit by 
the MCSM. (c) First-order extrapolation of the quadrupole 
moment of the 2f state. The results obtained by exact diago- 
nalization are also shown by the corresponding open symbols. 



t = 7 energy, —203.132 MeV. Nevertheless, the extrap- 
olated energy of the MCSM is -203.202 MeV, which is 
still much closer to the exact result than that of PHT. 
This good agreement provides us with a promising per- 
spective for its application to larger systems. 

Figure [3] shows the results of Of , 2f , Of , and Of states 
in order to discuss the behavior of the MCSM extrapola- 
tion concerning excited states and some observables. In 
Fig.[3ja), all of the MCSM-extrapolation results of these 
energies agree excellently with the exact ones in a unit 
of keV, too. Figures EJb) and (c) show the occupation 
numbers of the O/7 /2 orbit and the quadrupole moment of 
the 2f state by the MCSM and their first-order extrap- 
olations. In the case of these observables, a first order 
polynomial is appropriate for the extrapolation because 
the positive and negative contributions of the contam- 
ination of excited states cancel each other. Obviously, 
such cancellation does not occur in the case of energy 
eigenvalue. The first-order extrapolation for these ob- 
servables provides us with excellent improvement of the 
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FIG. 4: (Color online) Second-order extrapolations of the 
ground-state and 2 + energies of 64 Ge in the p/+gg/2-shell. 
The blue dashed line shows the first-order extrapolation of 
ground-state energy of PHT calculation. The notation is the 
same as that of Fig. [2] 



agreement with the exact value, while some other extrap- 
olation methods do not [l3j . 

Finally, we discuss the case of 64 Ge with pf + g 9 / 2 
model space in order to demonstrate the applicability of 
the present method to large-scale SM calculations. Its m- 
scheme dimension is 1.7 x 10 14 , which is roughly 10 3 times 
larger than the current limitation of the conventional di- 
agonalization method, ~ 10 11 . We adopt the PFG9B3 
effective interaction 18], which was used also in Ref . [l9j . 
In Fig. HI the result of the MCSM-extrapolation method 
shows stable behavior while the exact value is not avail- 
able. The 82 points for the ground state are obtained by 
the MCSM and fitted by a quadratic curve. The exci- 
tation energy of 2 + state is 0.95 MeV, which is close to 
the experimental value, 0.90 MeV [2(|. We also see the 
reasonable agreement between the ground-state energy of 



MCSM-extrapolation and that of the first-order extrapo- 
lation with PHT calculation. We point out that the PHT 
extrapolation is based on the four points (1 < t < 4), and 
the fitted line shows certain deviations from these points 
already, suggesting possible ambiguities. Note that the 
guess by the statistics of the nuclear level density is rather 
low, -306.7 MeV 0. 

In summary, we have proposed the MCSM- 
extrapolation method which provides us with accurate 
correction to the MCSM. Eq. ([5]) considerably reduces 
the computation time by orders of magnitudes to 
calculate the energy variance with deformed Slater 
determinants. The energy as a function of its variance 
is well fitted by a quadratic curve, and the result of the 
MCSM is improved down to a unit of keV especially in 
56 Ni case. We demonstrate that this method works quite 
well not only for energy eigenvalues, but also for other 
physical quantities of some low-lying states. By adopting 
the extrapolation method with the energy variance, we 
obtain a self-contained framework which removes the 
ambiguity of the energy convergence in the MCSM. We 
applied this framework also to large-scale shell-model 
problems, like the case of 64 Ge, which cannot be solved 
by existing conventional solvers. These results look 
quite promissing and encourages us to apply the present 
method to larger-scale problems. In such cases, the error 
estimation of the extrapolation method itself becomes 
important, and will be discussed in future publication. 
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